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In this general introduction to Ultra High Energy Cosmic Rays (UHECRs) we discuss the 
propagation of UHE protons and the GZK feature that is expected approaching 10 20 eV for 
homogeneously distributed sources. We also briefly present the effects of the propagation on 
other particles that can play the role of UHECRs. With the help of numerical simulations for 
the propgation of UHECRs, we show that the GZK feature cannot be accurately determined 
with the small sample of events with energies ~ 10 20 eV detected thus far by the largest two 
experiments, AGASA and HiRes. 

1 Introduction 

In the past ninety years of cosmic rays research there has been a constant search for the end of 
the cosmic-rays spectrum and it has long been thought that this end would be determined by the 
highest energy that cosmic accelerators might be able to achieve. Despite the continuous search, 
no end of the spectrum was found. In 1966, right after the discovery of the cosmic microwave 
background (CMB), it was understood^ that high energy protons would interact inelastically 
with the photons of the CMB and produce pions. For homogeneously distributed sources this 
would cause a flux suppression, called the GZK cutoff: for the first time the end of the cosmic 
ray spectrum was related to a physical process rather than to speculations on the nature of the 
accelerators. Moreover, for the first time, the end of the cosmic ray spectrum was predicted to 
be at a rather well defined energy, around 10 20 eV, where the so-called photo-pion production 
starts to be kinematically allowed. 

UHECRs can be of various nature and during their propagation over cosmological distances 
they suffer different kinds of energy losses. In this paper we consider most of the particles that 
could play the role of UHECRs and we review the processes affecting their propagation. In §2 we 
discuss the propagation of protons, heavy nuclei, photons and neutrinos. In §3 we show that the 
two largest experiments operating up to now, AGASA and HiRes, have a too small statistic to 



provide a conclusive answer about the presence or absence of the GZK feature in the UHECRs 
spectrum. We conclude in §4. 



2 Propagation of UHECRs in the cosmic photon background 

2. 1 Protons 

There are three sources of energy loss for ultra high energy protons propagating over cosmological 
distances: the expanding universe redshift, pair production (jyy — > pe + e~) and photo-pion 
production (jry — > ttN), each successively dominating as the proton energy increases. 

For protons the most important background is the CMB and the most important process 
is the photo-pion production in which a nucleon of sufficiently high energy sees, in its reference 
frame, the photons of the CMB blue-shifted to 7-rays above the threshold energy for photo-pion 
production, ^ ab - thr = m?r + m\j{2m N ) ~ 160 MeV. The cross-section for this process has a 
pronounced resonance just above threshold, corresponding to the production of an intermediate 
state A + that immediately decays into a nucleon and a pion, whereas in the limit of high energies 
it increases logarithmically with s = m 2 N + 2m^E^ h , giving rise to multiple pion production. 
For a background photon of energy e in the cosmic rest frame, defined as the frame in which the 
CMB is isotropic, the threshold energy £^ ab < thr translates into a corresponding threshold for the 
nucleon energy: 



^thr = r, ^ 2n rn N + ^ ) ~ 6.8-10 19 ( ^-^ ] ( —t— ) eV . (1) 
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Typical CMB photon energies are of the order of 10 -3 eV, giving a threshold value of a few tens 
of EeV (for a head-on collision). 

The interplay of this threshold with the Planck spectrum of the CMB photons produces 
a very steep, exponential, curve for the interaction length. The combination with the large 
inelasticity of the photo-pion interaction (the mean inelasticity goes from ~ 0.13 at threshold 
to ~ 0.5 at high energy, with large fluctuations) creates a very efficient and rapid mechanism 
to reduce the nucleon energy and makes the universe opaque to nucleons with energy above 
~10 20 eV on scales above ~100Mpc. The so-called GZK cutoff is due exactly to this: the flux 
at earth of nucleons with energy below threshold, say 5-10 19 eV, is due to contributions from 
(almost) all the universe, from Fig. ^ the loss length at this energy is of the order of 1 Gpc, 
while doubling the energy the loss length is reduced to 100 Mpc, and only a small portion of the 
universe contributes to the flux. Thus this change by a factor two in the energy changes the loss 
length by almost an order of magnitude, which translates in about the same ratio between the 
flux below 5-10 19 eV and above 10 20 eV if the sources have no luminosity evolution and no local 
over density and there is no magnetic field. 

Below ~ 6 • 10 19 eV the dominant loss mechanism for protons becomes the production of 
electron-positrons pairs on the CMB, p'j — » pe + e~ , down to the corresponding threshold: 



L , hi = ^£ (m7V + me ) ^4.8-10 17 f^-^j eV. (2) 



The interaction length for this process is much shorter than the one for pion production, but on 
the other hand the inelasticity is much lower, ~ 10~ 3 . This makes the pair production loss length 
of the order of Gpc (see Fig. The low inelasticity of pair production allows in calculations 
to treat this process as a continuous energy loss, whereas the pion production has to be treated 
as a discrete process due to its large inelasticity. 
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Figure 1: Solid line: loss length 
for photo-pion and photo-pair pro- 
duction for protons The 
dashed lines report the separate 
contributions of the two processes. 
The dotted line shows the redshift 
losses. 
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Figure 2: Modification factors as a function of the energy for many- 
source spectrum with 7 = 2.1 (solid lines) and 7 = 2.7 (dashed 
lines). The sources are uniformly distributed up to the indicated 
distances. After Ref.^1 



The last important mechanism which dominates near and below the pair production thresh- 
old is redshifting due to the expansion of the universe. Fig. ^ shows the loss lengths for pion 
and pair production as calculated in Ref.^. 

It is worth stressing that what has been named the GZK cutoff is in fact a feature® as 
the shape of the energy spectrum around 10 20 eV depends on many unknowns. The modifica- 
tions of the spectrum shape due to the above-mentioned loss processes was first investigated 
by Berezinsky and Grigorieva in Ref M They calculated the modification factor (basically the 
observed spectrum divided by the injection spectrum) for a uniform distribution of sources up to 
a maximum distance d max . Fig.|2]shows their results for sources without cosmological evolution, 
m = 0, for some values of the maximum distance of the sources. For large <i max , which is the 
case we are interested in, the spectrum shows a steepening followed by a flattening and then by 
a suppression. The flattening is due to the interplay between the features produced by the pair 
and pion production processes and it is an important feature for these spectra since it has a 
characteristic shape. There are claims that this feature has been observed in the experimental 
data^, although it is not yet clear if the feature in the data is due to this effect or if it is due to 
the transition between the galactic and extra-galactic components. 

It is important to stress what we said above: what is generically called GZK-cutoff is actually 
a feature as the spectrum does not end at 10 20 eV (see Fig. |2j), but has a flux suppression that 
depends on many details such as the injection spectrum of cosmic rays, the luminosity evolution 
of the sources, the local overdensity of sources and the magnetic field strength in the intergalactic 
medium. As an example, including the luminosity evolution makes the sources at high redshift 
brighter that the nearby ones and this enhances the flux suppression, while a local overdensity 
of sources has the opposite effect^ a flatter spectrum produces a lesser attenuation than a 
steeper one and the strength of the magnetic field in the intergalactic medium con produce 
many interesting features, see for example Ref.^. 



2.2 Heavy Nuclei 

For nuclei the situation is slightly different: the dominant loss process above about 10 19 eV is 
photodisintegration in the CMB and IR background (IRB) due to the giant dipole resonance, 
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Figure 3: Effective energy loss length for Fe photodisintegration off microwave (CMB), infrared (IR) and optical 
(OP) photons, as well as the total one (solid line) and the pair production loss length (PC). From Ref.^M 



except for light nuclei at very high energies^". Indeed, for a nucleus of mass number A and 
charge Z, the energy loss length for pion production is roughly the same one of a nucleon with 
identical Lorentz factor. This is due to the fact that the cross section for pion production is 
approximately proportional to the mass number A, while the inelasticity is proportional to 1/A. 
For pair production we got a different behavior because, while the inelasticity is proportional to 
1/A as before, the cross section is proportional to Z 2 resulting in an energy loss length lower by 
a factor A/Z 2 with respect to a proton with the same Lorentz factor. Since Z^A/2, the ratio 
of the photo-pair and photo-pion production increases roughly linearly with 

The cross sections for photodisintegration <7A,i(e') contains essentially two regimes depending 
on e', the photon energy in the nucleus rest frame. At e' < 30MeV there is the domain of the 
giant dipole resonance and the disintegration proceeds mainly by the emission of one or two 
nucleons. At higher energies, the cross section is dominated by multi- nucleon emission for heavy 
nuclei and is approximately flat up to e'~150MeV. A useful quantity to estimate the energy 
loss rate by photodisintegration is given by the effective rate: 



For photodisintegration, the average fractional energy loss results equal to the fractional loss in 
mass number of the nucleus, E~ 1 dE/dt = A~ 1 dA/dt, because the nucleon emission is isotropic 
in the rest frame of the nucleus. Therefore during the photodisintegration process the Lorentz 
factor of the nucleus is conserved, unlike the cases of pair and pion production which involve the 
creation of new particles that carry away energy. The energy loss time for photodisintegration is 
then A/R e ^. Fig. |31 shows separately the different contributions to this quantity from CMB, IR 
and optical photons for Fe nuclei, together with the total one (solid line) and the pair creation 
loss length. 

It is apparent that the optical background has no relevant effect, that the IR one dominates 
the photodisintegration processes below 10 20 eV and the CMB dominates above 10 20 eV. The 
pair creation rate is relevant for Fe energies 4-10 19 eV -r- 2-10 20 eV (7 factors ~ (1 4- 4) • 10 9 ), 
for which the typical CMB photon energy in the rest frame of the nucleus is above threshold 
(> 1 MeV) but still well below the peak of the giant resonance (~ 10 4- 20MeV). The effect of 
pair creation losses is to reduce the 7 factor of the nucleus, obviously leaving A unchanged^. 
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Figure 4: Average energy as a function of the propagation dis- 
tance for particles that started as Iron nuclei with the indicated 
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Lorentz factors. From Ref.^- 1 . 



Figure 5: Interaction lengths (dashed lines) 
and energy attenuation lengths (solid lines) 
of 7-rays in the CMB (thin lines) and in the 
combined CMB and URB (thick lines). The 
interactions taken into account are single and 



double pair production. From Ref. 



We should not get fooled by the loss lengths in Fig. |31 into thinking that if the loss length 
for a Fe nucleus of 10 20 eV is 500 Mpc, then we can receive on Earth a Fe nucleus that started 
many hundred Mpc away. This is because the corresponding interaction length is more than an 
order of magnitude shorter and after every interaction the nucleus becomes lighter and lighter 
and alongwith this the loss length for photodisintegration becomes shorter and shorter. The 
net result'-", as can be seen in Fig. 0J is that after 10 Mpc all the energies are below 2-10 20 eV 
and after 100 Mpc they are below 10 20 eV. 



2.3 Photons 

As in the case of UHE nucleons and nuclei, the propagation of UHE photons (and electrons/pos- 
itrons) is also governed by their interaction with the cosmic photon background. The dominant 
interaction processes in this case are the attenuation of UHE photons due to pair production 
(PP) on the background photons (77^ — > e + e~), and inverse Compton scattering (ICS) of the 
electrons (positrons) on the background photons. 

The 7-ray threshold energy for PP on a background photon of energy e is 

S thr = ^^2.6-10 11 ^)~ 1 eV, (4) 

whereas ICS has no threshold. In the high-energy limit, the total cross sections for PP and ICS 
are: 

3 

app ~ 2cr IC s ^ ^o- T (m 2 e /s) ln(s/2m 2 ) (s > m 2 ) . (5) 

For s <C m 2 , cries approaches the Thomson cross section <tt = 8ira 2 /3m 2 (a is the fine structure 
constant), whereas crpp peaks near the threshold. Therefore, the most efficient targets for 
electrons and 7-rays of energy E are background photons of energy e ~ m 2 /E. For UHE this 
corresponds to e ~< 10~ 6 eV ~ 100 MHz. Thus, radio background photons play an important 
role in UHE 7-ray propagation through extragalactic spaced Unfortunately, the universal radio 
background (URB) is not very well known, mostly because it is difficult to disentangle the 
Galactic and extragalactic components. 



In the extreme Klein-Nishina limit, s ^> m 2 , either the electron or the positron produced 
in the process 77^ — > e + e~ carries most of the energy of the initial UHE photon. This leading 
electron can then undergo ICS whose inelasticity (relative to the electron) is close to 1 in the 
Klein-Nishina limit. As a consequence, the upscattered photon which is now the leading particle 
after this two-step cycle still carries most of the energy of the original 7-ray, and can initiate 
a fresh cycle of PP and ICS interactions. This leads to the development of an electromagnetic 
(EM) cascade which plays an important role in the resulting observable 7-ray spectra. An 
important consequence of the EM cascade development is that the effective penetration depth 
of the EM cascade, which can be characterized by the energy attenuation length of the leading 
particle (photon or electron/positron), is considerably greater than just the interaction length 
(see Fig. Hjfl. 

EM cascades play an important role particularly in some exotic models of UHECR origin 
such as collapse or annihilation of topological defects in which the UHECR injection spectrum 
is predicted to be dominated by 7-rays. But, even if only UHE nucleons and nuclei are produced 
in the first place, for example via conventional shock acceleration, EM cascades can be produced 
by the secondaries coming from the decay of pi ons which are created in interactions of UHE 
nucleons with the low energy photon background!*^*. 

Most of the energy of fully developed EM cascades ends up below ~ 100 GeV where it is 
constrained by measurements of the diffuse 7-ray flux. Flux predictions involving EM cascades 
are therefore an important source of constraints of UHE energy injection on cosmological scales. 

It should be mentioned that the development of EM cascades depends sensitively on the 
strength of the extragalactic magnetic fields (EGMFs) which is rather uncertain. The EGMF 
typically inhibits the cascade development because of the synchrotron cooling of the e + e~ pairs 
produced in the PP process. The energy lost through synchrotron radiation does not, however, 
disappear; rather, it reappears at lower energies and can even initiate fresh EM cascades. 

2.4 Neutrinos 

The propagation of UHE neutrinos is governed mainly by their interactions with the relic neu- 
trino background (RNB). The interaction energies are typically smaller than electroweak energies 
even for UHE neutrinos and then the cross sections are given by the Standard Model of elec- 
troweak interactions which are well confirmed experimentally. Physics beyond the Standard 
Model is not expected to play a significant role in UHE neutrino interactions with the low- 
energy relic backgrounds. Despite the neutrino-neutrino cross section are at least a few order 
of magnitude smaller than the neutrino-nucleon ones, the latter interactions are negligible com- 
pared to interactions with the RNB because the RNB particle density, ~ 100 cm -3 per family, 
is about 10 orders of magnitude larger than the baryon density. 

The vV annihilation mean free path is of the order of \ v = [n v a v -^)~ x ~ 4-10 28 cm, just above 
the present size of the horizon (Hq 1 ~ 10 28 cm). The neutrino is the only known stable particle 
that can propagate through the universe essentially uninhibited even at the highest energies. 
This has lead to the speculation that neutrinos could be indeed the super-GZK primaries. 
However, in the Standard Model a neutrino incident vertically in the atmosphere would pass 
through it uninhibited, never initiating an extensive air shower. Consequently, for these scenarios 
to work, one has to postulate new interactions so that these neutrinos acquire a strong cross 
section above 10 20 eV. 

An interesting situation arises if the RNB consists of massive neutrinos with m u ~ 1 eV: 
such neutrinos would constitute hot dark matter which is expected to cluster, for example, in 
galaxy clusters. This would potentially increase the interaction probability for any neutrino 
of energy within the width of the Z resonance at E = M 2 z /2m v = 4 • 10 21 (eV/m^) eV. It 
has been suggested that the stable end products of the Z-bursts induced at close-by distances 



(~< 50Mpc) from Earth may explain the highest energy cosmic rays^*^. The problem with these 
proposals is however that they require a very high flux of UHE neutrinos to begin with and this 
makes Z-burst above GZK energies more likely to play a role in the context of non-accelerating 
scenarios. For further information see Ref.^and references therein. 

It is important to point out that the only conventional/assured source of UHE neutrinos is 
the GZK effect itself. The neutrinos are the result of the decay of the pions produced in the p"f 
interaction. The flux howe ver is not very high and the detection is quite difficult. For further 
informations see Refs.^LL^l 

3 AGASA and HiRes: is there a discrepancy? 

AGASA and HiRes are, up to now, the two experiments with the larger exposure for the detection 
of UHECRs. They reported however apparently conflicting results. The two reported spectra 
appear: 1) to have a systematic offset at low energy and 2) to differ above 10 20 eV where AGASA 
shown no hint of the GZK-suppression whereas HiRes seems to be consistent with it. It has been 
shown^l that a systematic overestimate of the AGASA energies by 15% and a corresponding 
underestimate of the HiRes energies by the same amount woul d in fact bring the two data sets 
in a much better ag reement in the region below 10 20 eV. In Ref.^^we applied our Monte Carlo 
simulation! I I I to investigate the discrepancies at high energy and we found that: 

• assuming a uniform distribution of sources, the AGASA spectrum is reproduced, in a con- 
ventional scenario where the average spectrum has the GZK suppression, with a probability 
of - 6-KT 4 (~ 3a). 

• assuming the presence of the 15% systematic error, the shifted AGASA data are reproduced 
with a probability of ~ 6-KT 3 4- 1CT 2 (2.5a 4- 2a). 

• the HiRes data are reproduced, in a scenario without a GZK suppression 51 , with a proba- 
bility of ~ 2%, (~ 2a). 

It is important to stress that in order to properly do these calculations one has to take into 
account the statistical error in the energy determination. Due to the steeply falling spectrum and 
to the expected change of slope around 10 20 eV the statistical error in the energy determination 
changes the average expected number of events above 10 2 ~eV by ~ 1 and with the present 
limited statistics even a difference of one event is important^! 

In Fig. El we plot the spectra of some of the simulated AGASA realizations that produced 
11 or more events above 10 20 eV. It is striking the resemblance of the simulated spectra to 
the AGASA one: all of them show no evidence of the GZK suppression. This shows that the 
AGASA spectrum is far from being impossible, even if the average cosmic ray spectrum can be 
expected to show a GZK feature. 

From the above points we can conclude that neither AGASA nor HiRes have enough statis- 
tical power to prove the presence or absence of the GZK feature in the spectrum of UHECRs. A 
new generation of experiments is needed to finally provide a conclusive answer to this question. 

4 Conclusions 

We considered several particles that can possibly play the role of UHECRs and we discussed 
the energy loss processes that affect their propagation over cosmological distances. We showed 
that for all considered particles the energy losses above ~ 10 20 eV become so severe that they 
cannot propagate over distances larger than about 100 Mpc without reducing their energy below 



a To mimic a scenario without GZK suppression we used the AGASA dataset as template. 
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Figure 6: In the above panels we plot 4 of the 18 simulations that have 11 or more events above 10 20 eV. The 
black crosses are the simulation results. The red dots with errorbars are the AGASA data superimposed for 
comparison. The errorbars in the AGASA data are slightly shifted left to avoid covering up the black crosses. 

10 20 eV. The effect of these energy losses on the spectrum of UHECRs is the so called GZK- 
suppression, due to the fact that below 10 20 eV almost all the universe is contributing to the 
observed flux whereas above 10 20 eV we receive contributions only from sources not too far away, 
~ 100 Mpc. We stressed the point that this is not a cutoff, but a feature since the spectrum 
does not end at 10 20 eV, but it is only suppressed, and the amount of this suppression depends 
on many unknowns such as: the luminosity evolution of the sources, their local overdensity and 
the magnetic field strength in the intergalactic medium. 

We showed that the present sets of data are not enough to determine whether this GZK 
suppression is present or not in the observed spectrum and that we need a new generation of 
experiments to have a conclusive answer to this question. 
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